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ABSTRACT 


The  generation  of  surface  waves  in  a  fluid  caused  by  a 
partially  submerged  body  is  modeled  mathematically  as  a  boundary 
value  problem  for  the  Laplacian  with  Neumann  data  on  the  bottom 
of  the  fluid  container  and  on  the  body,  a  linearized  free  sur¬ 
face  condition  and  a  radiation  condition  at  infinity.  Fritz 
John  (Comm.  Pure  Appl.  Math. ,  III,  1950)  shoved  that  the  problem 
had  a  unique  solution  only  under  rather  reatrictive  assumptions 
on  the  body  geometry  and  reformulated  the  problem  as  an  integral 
equation  which  was  not  uniquely  solvable  at  certain  irregular 
frequencies.  In  the  present  work,  uniqueness  is  established  for 
more  general  geometries,  allowing  corners  and  nonnormal  inter¬ 
sections  of  the  body  with  the  free  surface.  Also  presented  are 
two  methods  of  modifying  the  integral  equation  so  that  it  is 
uniquely  solvable  for  all  frequencies.  One  method  involves 
introducing  an  additional  integral  on  the  waterplane,  while  in 
the  second  method,  an  additional  integral  term  is  added  to  the 
equation  which  remains  an  equation  only  over  the  submerged  por¬ 
tion  of  the  body. 
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INTRODUCTION 

It  has  been  thirty  years  since  Fritz  John  published  a  pair  of  papers  with  a 
similar  title1*’2  which  still  underlie  present  thinking.  In  Reference  2,  which  must 
be  regarded  as  a  tour  de  force  of  classical  applied  mathematics,  John  analyzed  simple 
harmonic  motion  of  the  fluid  in  which  an  impenetrable  body  is  partially  immersed. 

He  formulated  the  problem  mathematically  as  a  boundary  value  problem  for  Laplace's 
equation  with  appropriate  boundary  conditions  on  the  body,  the  free  surface,  the 
bottom  of  the  fluid  container,  and  the  radiation  condition.  Under  certain  restric¬ 
tions  on  the  body  shape,  he  proved  the  existence  of  a  unique  solution  and  formulated 

*A  complete  listing  of  references  is  given  on  page  57. 
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an  integral  equation  for  the  velocity  potential  evaluated  on  the  body  surface.  This 

equation  is  not  uniquely  solvable  at  a  set  of  irregular  frequencies  as  John  pointed 

out  and  this  has  plagued  researchers  in  this  area  ever  since. 

John  based  his  work  on  classical  potential  theory  even  though  the  physical 

domain  had  a  corner  at  the  intersection  of  the  body  and  the  free  surface.  Thus,  he 

looked  for  classical  solutions,  continuously  differentiable  up  to  and  including  the 

boundary,  although  he  did  not  show  that  the  solution,  the  existence  of  which  he 

proved,  actually  satisfied  this  property.  This  is  not  surprising  in  the  light  of 
3 

subsequent  work  which  shows,  in  fact,  that  such  smooth  solutions  do  not  exist.  How¬ 
ever,  by  employing  the  fundamental  work  on  potential  theory  for  irregular  domains, 
developed  by  Burago  et  al..  Krai,  and  Uendland,  it  is  possible  to  formulate  the 
problem  in  a  more  appropriate  mathematical  setting  and  to  show  the  existence  of  a 
unique  solution  under  less  restrictive  assumptions  on  the  obstacle  shape.  This  is 
done  in  the  present  report  where  uniqueness  and  existence  are  established  for  piece- 
wise  smooth  bodies.  If  we  permit  the  wave  number  to  have  a  nonzero  imaginary  part, 
no  further  restrictions  are  required,  whereas  we  do  restrict  the  bodies  to  those 
whose  projection  onto  the  free  surface  coincides  with  the  waterplane  area  if  the  wave 
number  has  real  values.  In  either  case,  we  relax  the  conditions  John  imposed  and 
permit  nonnormal  intersections  of  the  body  with  the  free  surface  as  well  as  bodies 
with  corners. 

We  also  present  two  methods  of  resolving  the  problem  of  irregular  frequencies. 
One  is  essentially  a  simplified  version  of  the  method  proposed  by  Chang  and  Pien^  in 
which  the  obstacle  surface  is  augmented  by  the  waterplane  to  obtain  an  integral  equa¬ 
tion  on  a  closed  surface  soluble  at  all  frequencies.  The  second  method  is  modeled 

after  a  treatment  of  exterior  problems  for  the  Helmholtz  equation  where  similar  prob- 

8  9  10 

lems  occur;  Burton  and  Miller,  Kleinman  and  Roach,  and  Angell  and  Kleinman. 

Here  the  integral  operator  is  augmented  rather  than  its  domain  and  range  which 
remain  the  submerged  portion  of  the  body.  An  alternate  method  for  removing  the  ir¬ 
regular  frequencies  has  been  given  by  Uracil,**  using  a  modified  Green's  function  as 
12 

proposed  by  Jones,  to  handle  a  similar  problem  in  acoustic  scattering. 

NOTATION  AND  FORMULATION 

As  illustrated  in  Figures  1  and  2,  the  geometry  of  our  problem  is  described  in 
the  following  manner.  Let  D_  be  a  connected  bounded  open  set  in  three-dimensional 

3  * 

space,  ft  ,  whose  boundary  consists  of  C^,  a  closed  bounded  connected  region  in  the 
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Figure  1  -  Floating  Body  -  Three-Dimensional  View 


Figure  2  -  Floating  Body  -  Section  in  (p,y)  Plane 


x-z  plane  and  Cq,  a  piecewise  smooth  surface  in  the  lower  half  space,  y  <  0,  in  the 

following  sense.  The  boundary  3D_  *  t0  b®  piecewise  Liapunov,  that  is,  con- 

posed  of  a  finite  number  of  segments  each  of  which  lies  on  a  Liapunov  surface  (on 
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which  there  is  a  Holder  continuous  normal);  see,  e.g.,  Gunter  for  a  precise  defini¬ 
tion  of  Liapunov  surfaces.  Moreover,  3D_  satisfies  a  two-sided  cone  condition;  i.e., 
there  exist  a,  h  >  0,  such  that  each  point  of  3D_  may  be  taken  as  the  vertex  of  two 
cones  of  height  h  and  vertex  angle  o,  one  lying  entirely  in  D_  and  one  in  the  com¬ 
plement  of  the  closure  of  D_.  A  rectangular  coordinate  system  is  oriented  with  the 
origin  in  and  D_  lying  below  the  y  *  0  plane.  Let  C^,  the  free  surface,  be  the 
closure  of  the  complement  of  in  the  y  ■  0  plane,  and  let  Cg  be  the  plane  y  *  -h. 

Denote  by  DL  the  unbounded  domain  bounded  by  C,,  C  ,  and  C_.  Also  let  (f  be  the  re- 
flection  of  Cq  in  the  x-z  plane  and  D  be  the  reflection  of  D_.  Unless  otherwise 

<V 

noted,  the  normal  n  points  into  D,  from  C  ,  Cc,  C_,  and  into  D  from  C  .  We  denote 

r  +  o  f  B  w 

points  (vectors)  by  p  •  (xp,yp,zp)  and  q  *  (x^.y^.z^)  with  cylindrical  coordinates 
p  *  (Pp,0p,yp)  a"d  the  subscripts  will  be  omitted  if  there  is  no  danger  of  confusion. 

We  formulate  the  floating  body  problem  as  follows: 

Find  a  function 

<J>(p)  e  C2(D+)nC°(D+u3D+)  (la) 

where  Cm  (ft)  is  the  space  of  functions  with  m  continuous  derivatives  on  ft,  such  that, 

V2$  -  0  p  e  D+  (lb) 


5  A 

35* 


w  -  + 


k* 


&P  m  W  m  n 
35  37  0 


d<t> 

35 


V(p) 


p  e  Cf  (lc) 

P  e  CB  (Id) 

p  e  CQ  (le) 
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where  V(p)  is  integrable  on  C  , 

o 


where  kQ  is  the  root  with  largest  real  part  of  the  equation 

kQ  sinh  kQh  -  k  cosh  kQh  (2) 

Some  properties  of  kn  are  listed  in  Appendix  A;  here  we  only  remark  that  there  always 
is  such  a  kQ  with  a  real  part  that  is  strictly  positive. 

Finally  we  add  a  condition  on  3$/3n  which  will  enable  us  to  employ  Green's 
theorem,  namely,  that  3<J>/3n  defines  a  distribution  in  C°*(3D+)  (the  space  of  func¬ 
tions  of  bounded  variation)  in  the  sense  of  boundary  flow  (see  Burago,  Haz'ja  and 
Sapoznlkhova,  Wendland,6  and  Klelnman  and  Wendland1*  for  a  more  detailed  discussion 
of  this  concept).  Roughly  this  means  that 

lim  f  u  ■  ds 

3D  -*-3D.  I  3n 

m  +  J 

3D 

m 

defines  a  bounded  linear  functional  on  C° (3D, )  (continuous  functions  on  3D.)  where 
«»  3  +  +  , 

u  £  CQ  (  J<  )  (infinitely  differentiable  functions  with  compact  support  in  X)  and 

3Dm  is  a  family  of  smooth  surfaces  in  D+  converging  to  3D+.  We  note  that,  if  (3$/3n) 

£  Ljoc  (3D+) ,  it  will  have  boundary  flow.  Alternatively,  we  could  require  that 

♦  (p)  e  H^oc  (D+),  that  is,  J  { |<f  j  |  ^}dt  <  *  for  every  bounded  subdomain  ft  c  D+. 

This  would  also  enable  us  to  employ  Green's  theorem.  Klelnman  and  Wendland14  have 
shown  that,  if  <p  has  boundary  flow,  then  it  is  in  H1qc (D+) .  Hie  given  data  V  will 
be  interpreted  as  defining  a  boundary  flow  which  it  will  do  if,  for  example,  V  e 
1>  <CQ) .  Here,  L  (ft)  is  the  space  of  (Le  be  ague)  Integrable  functions  on  ft;  H'L(ft)  la 
the  space  of  functions  which,  together  with  their  first  derivatives,  are  in  Ll(ft); 
and  l|oc(Q)  is  the  space  of  functions  integrable  over  every  bounded  subdomain  of  ft. 


UNIQUENESS 

The  first  important  property  we  establish  concerns  conditions  under  which  the 

problem.  Equations  (la)- (If),  has  a  unique  solution.  We  remark  that  John  proved 

2 

uniqueness  by  requiring  that  Cq  was  smooth  (C  )  and  that  rays  perpendicular  from  Cfi 

intersect  C  at  most  once.  In  fact,  he  requires  C  (or  C  C  )  to  be  convex.  Here 
o  o  ou  o 

we  will  establish  two  versions  of  the  uniqueness  proof.  In  one  we  require  the  waves 

to  be  dispersive,  that  is,  Im  k  >  0,  in  which  case  uniqueness  is  established  for  the 

piecewise  smooth  boundaries  described  earlier.  The  second  deals  with  nondispersive 

waves,  Im  k  *  0,  and  requires  an  additional  assumption  that  vertical  rays  from 

(not  Cg  as  in  John's  proof)  do  not  contain  points  in  D_.  In  either  case,  we  make  no 

convexity  requirements  and  permit  nonnormal  intersections  of  C  and  C,.  We  note  that 

or  ^5 

the 'assumption  of  a  small  imaginary  part  for  k  has  ample  precedent  (e.g.,  Llghthill 

and  Noblesse. ^  Also,  it  should  be  noted  that  Lenoir  and  Martin*-2  have  presented  a 

uniqueness  proof  for  real  k  in  the  case  of  Infinite  depth  which  requires  no  geometric 

restrictions.  Their  method,  however,  does  not  appear  to  be  immediately  extendable 

18 

to  finite  depth.  In  fact,  Ursell  has  shown  that  their  proof  is  incorrect. 
Uniqueness  Theorem: 

If  Cq  is  piecewise  smooth  and  Im  k  >  0,  the  homogeneous  floating  body  problem 
has  only  the  trlval  solution;  that  is,  if 


V2<j>  -  0,  P  e  D+ 

3$  _  n  _  _  r 

3n  °»  P  E  o 

~  ■  0,  p  e 

-  lk0*  ■  °  (ri7i) 

(3) 

!y  "  °*  P  E  °B 

♦  e  C2  (D+)n  C°  (D+U3D+) 

and  3$/3n  exists  as  a  boundary  flow  even  though  it  is  not  defined  at  corners,  then 
4>  =  0.  The  theorem  remains  true  if  Im  k  ■  0,  and  vertical  rays  from  do  not 
Intersect  D  . 
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Proof : 

'  Assume  4>  satisfies  Equation  (3).  Since  <t>  is  continuous  in  C°  (D+(J3D+)  and  has 
a  normal  derivative  in  the  sense  of  boundary  flow,  we  may  still  apply  the  divergence 
theorem  (Burago  and  Maz'ja)19  to  <j>V$  -  $V<J>  in  D+u  DR,  where  DR  -  {pjp<R}  is  a 
cylinder  of  radius  R;  and  $  denotes  the  complex  conjugate,  thus  obtaining 


<*v2?v*>dT .  -j  {♦!!-?  ii} 


°+nDR 


Cf  UCoUCBU3DR 


where  h  points  into  D+.  Because  <}>  satisfies  Laplace's  equation,  the  integral  on  the 
left  vanishes  and  the  boundary  conditions  cause  the  Integrals  over  Cq  and  CR  to 
vanish  as  well.  Employing  the  boundary  condition  on  and  the  fact  that  on  3DR, 
3/3n  »  -  3/3p,  we  obtain 


(E-w  |*i2  d.  +  J 


The  radiation  condition  implies  that  on  3DR 


$  -  ‘V  2  •  (t  -  “„♦)  (f +  *M) 


+  lyi2  + e,  (♦  H + ?  H) + (*  H  -  ♦  H)  <*> 


hence. 


•  <>(£)■  °(i) 


H ' 5  Ip  ■  iei;  { Ip  +  lfc0*t2+I"  ko  (♦  If  +  *  If)}  +  ° (I) 
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m4 


Substituting  in  Equation  (5)  we  obtain 


f  (t"k>  '*'2  d® +  (|lp  +  iv'2+im  k0  (♦  H + *  !£)}d8 " o(1) 

C.  °3D.  1  ’ 


as  R  -*•  00  (8) 


or 


2  Ink 


J  I*'2  ds  +  R^J  {|lp  +  iVl2+Iako  (*  If  ^  H)}  ds  “  o(1) 

Cf  3DR  as  R 


as  R  -*•  00  (9) 


From  property  4  in  Appendix  A  we  know  that  if  k  j*  0,  then  RekQ  >  0.  Equation  (9) 
may  be  written  as 

2  Im  k  I"  |<j>|2  ds  J  j(Reko)2  |<J>|2  +  (im  kQ  R&p  +  Re 

r  °  an  ' 


3D 


+  (im  kQ  Im(J>+Im  ||) 


ds  »  o(l) 


(10) 


Since  Im  k  >  0,  then  all  the  terms  on  the  left  are  nonnegative,  and  it  follows  that 
individually  they  vanish  as  R  00 .  Thus, 


J  |<fr|2  ds  -  o(l)  =>  ♦  -  o  (^72) 

9DR 


(11) 


which,  together  with  the  remaining  relations,  ensures  that 


Ip "  °  (rb?) 


(12) 


and 


8 


I 


ds  -  o(l) 


03) 


But  the  last  expression  Is  Independent  of  R,  hence 


In  k 


ds  ■  0 


(14) 


Equations  (11)  and  (12)  Imply  that  $  3<fr/3n  -*■  0  on  3D^  as  R  hence  we  nay  apply 
the  divergence  theorem  to  in  D+  ,  obtaining 


-/  , 

CoUCfUCB 


|V4»|2  dT 


(15) 


which,  with  the  boundary  conditions,  reduces  to 


k 


dx 


(16) 


If  In  k  >  0,  then  Equation  (14)  implies  that 


ds  »  0  and. 


therefore. 


J  1 76 1 2  dx  -  0 

Thus  6  must  be  constant  and,  by  Equation  (14),  the  constant  nust  be  aero  and  the 
theorem  is  proven.  If  In  k  -  0,  Equation  (14)  gives  no  information  and  Equation 
(16),  although  still  valid  does  not  imply  that  ♦■O  since  k  has  the  "wrong"  sign. 
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4 


To  prove  uniqueness,  when  Is  k  ■  0  m  follow  John's  argument  and  lapose  the 
additional  restriction  that  every  vertical  ray  fron  the  free  surface  contains  no 
points  of  D_;  that  is,  if  (x,0,z)  e  Cf  then  (x,y,z)  t  D_,  for  -h  <  y  <  0.  This  re¬ 
striction  means  that  the  projection  of  D_  on  the  y  ■  0  plane  is  Cw,  the  waterplane, 
as  illustrated  in  Figure  3a.  We  introduce  additional 


Cw  C, 


Figure  3a  -  Permissible  Configuration  Figure  3b  -  Impermissible  Configuration 

Figure  3  -  Restricted  Geometry  for  Im  k  -  0 

notation  setting  Dj  to  be  that  portion  of  D+  that  may  be  reached  by  vertical  rays 
from  Cf  and,  following  John,  we  have  the  representation 


-2T  an(x,z)  cosh  kn(y+h)  for  p  e  (17) 

n-0 


where  k  are  the  roots  of  Equation  (2)  with  nonnegative  imaginary  part.  Using  Equa- 
n 

tions  (17)  and  (2)  we  may  show,  by  elementary  Integration,  that 


I  h  /  •inh  2koh\ 

J  4>(x,y,z)  cosh  kQ(y+h)  dy  -  j  •„<*.*>  ^1+  2k^h ) 

•h 


(18) 


With  the  Schwarz  inequality  we  find 


I 


-  lim  2nR  ^1/  *H|nl  (UoR) ' 


from  which  we  infer  that 


3n-° 


for'  tevery  n  (25) 


hence. 


ao  "  ° 


$  cosh  k  (y+h>  dy  -  0 


Integration  by  parts  leads  to 


hence. 


f° 

t.O.z)  sinh  kQh  -J  4>y  sinh  kQ(y+h)  dy 


|*(..0..>  .lot.  *0M2  if  I*,!2  *>[  •l"h2  ko<y+h)  * 


Since,  froa  Equation  (2), 


L 


ainn  ko(yfh)  dy  - 


■lnh  k  h  coeh  k  h 
O  P 

2k 


h 

2 


•infa  k  h 
o 

2k 


h 

2 


(29) 


we  have,  with  Equation  (28), 


2k  !♦(*,<>,«) j 2 


*T 


Integrating  over  we  obtain 


2k 


|9|2  da 


1 99 1 2  dx 


and,  with  Equation  (16), 


or 


|79|2  dr  <  0 


(») 


(31) 


(32) 


(33) 


froa  which  we  Infer,  alnce  both  teraa  on  the  left  are  nonnegative,  that  V9  -  0  in  D_, 
hence  9  -  constant .  But,  then  39/bn  -  0  and  the  free  surface  condition  tallies  that 
9*0  coapleting  the  proof.  We  note  that  the  proof  would  be  essentially  unchanged 
if,  inatead  of  39/dn  «  0,  we  had  9  a  0  on  CQ. 
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INTEGRAL  EQUATION  FORMULATION  ... 

Again  we  follow  John  and  introduce  the  Green's  function  for  the  limiting  case 
when  D_  is  empty  and  C{  is  the  entire  y  -  0  plana.  Thus,  define 


Y(p,q) 


Z 


r  <ln-k2>Hi1> <k„T> 


s  .  2  2 

hk  -bk  4k 
n-0  n 


cosh  k  (y  +h)  cosh  k  (y  +h) 
n  p  n  q 


00 


J  (MR)  cosh  y(y<  +  h)  H  HP  ±  *  Binh  py> 

o'^  J  y  sinh  yh  -  k  cosh  yh 


where  y>  *  max  {yp,y^}  £  0  and  y<  ■  min  {y^.y^}  i  -h; 


R  »  /(V*q)2+<V*q>2  (35) 

Here  kR  are  the  roots  of  Equation  (2)  with  nonnegative  real  and  imaginary  parts, 

and  the  contour  passes  below  the  zeros,  y  -  k  ,  of  the  denominator.  It  is  not  diffl- 

n 

cult  to  see  that  y  satisfies  the  radiation  condition  and  boundary  conditions. 


|l-0,  y--h  and  |^-ky*0,  y-0 


John  has  also  shown  that  the  singular  behavior  of  y  is  given  by 


y - - — r - - - 7 - +  H(p,q) 

27r|p-q|  2n  J  p-q  |  2tt  j  p-<Jx 


where  H  »  0(An|p-qrj)  0  (in  |  p-q.J  ) 

T  "  <WV 

qi  ’  (xq*“V2h,*q) 


He  now  state  Green's  theorem  for  solutions  for  the  floating  body  problem  using  the 
Green's  function  defined  in  Equation  (34). 


If 


V  $  -  0,  P  e  D+ 


fc  -  -  •  (S) 


|y  "  **  *  0*  P  E  Cf 
In  k  >  0 

and  3<f>/3y  exists  at  least  in  the  sense  of  boundary  flow,  then 


I 


■  Y(p»q)  Ur 

q  q 


ds  “  -  a(p)  <t>  (p) 


(38) 


where  n  points  into  D+  and 


f  ^L- 

J  3n 


a(p)  :■  lim 

£"°  [3B£(P)]nD+  q 


ds 


2*  p  e  D+UCfU  CB 

1,  p  e  smooth  points  of  C 


(39) 


0,  p  e  D 


-UCw 


The  fact  that  a(p)  -  2  on  and  rather  than  1  is  a  consequence  of  the  form  of  y. 

Equation  (37).  Here  B£(p)  denotes  a  sphere  of  radius  c  with  center  at  p,  the  normal 

points  away  from  B£ (p)  and  C'  and  are  primed  to  denote  that  the  boundary  points 

between  C,  and  C  are  not  included.  It  should  be  noted  that  a(p)  is  defined  at  non- 
t  w 

smooth  as  well  as  smooth  points  and  is  a  measure  of  a  normalised  solid  angle,  al¬ 
though  the  explicit  evaluation  in  Equation  (39)  holds  only  at  smooth  points. 

We  also  state  a  version  of  Green's  theorem  for  solutions  of  Laplace's  equation 
in  D_:  If  V  -  0,  p  e  D_  and  30/3n  exists  at  least  as  an  interior  boundary  flow 
(approximating  surfaces  lie  in  D_),  then 
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c  c 

ou  w 


T(al  2Y.(Pi<1>.  .  v  3L 


-  [2-a(p)J  <t> (p> 


where  points  away  from  D  . 

Since  3y/3n  =  ky  on  C  ,  this  may  be  written 

H  * 

j  I**’’  ^  >  M^j  V(P,q)  -O’)  d., 


(2-a(p)  ]  <Hp) 


In  particular,  if  $  *  1,  we  get  a  form  of  Gauss'  Integral 


J  3„q  , 


*‘f 


2  -  a(p) 


It  should  be  kept  in  mind  that  both  Equations  (38)  and  (42)  hold  for  points  p  in  the 
slab  £  yp  £  0  since,  as  yet,  we  have  not  defined  y  outside  of  this  region.  With 
the  boundary  condition.  Equation  (le) ,  and  Equation  (38),  we  have  the  boundary 
integral  equation 


a(P)  <1><P)  +  J  <Kq)  ^><L'  da  -  f  V(q)  y(p,q) 
_  a  ^  J 


ds  for  p  £  C  (43) 
q  r  o 


where  a(p)  =  1  at  all  smooth  points  of  C  .  Alternatively,  we  could  multiply  Equa¬ 
tion  (42)  by  <J>(p)  and  add  it  to  Equation  (38),  and  use  Equation  (le),  to  obtain 

2  <J>(p)  +  J  ['Kq)-<Kp)]  dsq  ~  k$(p)  J  Y(p,q)  dsq 


I  V(q)  Y(p.q)  dg 

v  q 

C 

o 


for  p  t  D+u  3  D+  (44) 


* 


This  is  a  form  derived  by  Noblesse16  and  is  the  analog  of  the  formulation  for  the 
Helmholtz  equation  derived  by  Ahner  and  Kleinman.  Observe  that  the  representation 

does  not  change  form  when  p  is  on  the  boundary  of  even  if  p  is  a  corner  point. 

If  we  restrict  p  to  lie  on  Cq,  both  Equations  (43)  and  (44)  are  integral  equations 
for  the  unknown  values  of  ♦  on  Co<  As  John  showed,  there  exist  irregular  values  of 
k,  that  is,  those  values  for  which 

o(p)  u(p)  +  J  u(q)  ds  ■  0  for  p  e  C  (45) 

JC  *  ° 

o 

has  nontrivial  solutions.  Ve  note  that  if  u  is  a  solution  of  Equation  (45),  it  is 
also  a  solution  of  the  homogeneous  form  of  Equation  (44).  Note  also  that  since  a(p) 
■  1  almost  everywhere  on  Cq,  Equation  (45)  may  also  be  written 


a(p)  u(p)  +  f  a(q)  u(q)  *  0  for  P  e  CD  (46) 

Co 

and  so,  even  if  CQ  is  only  piecewise  smooth,  the  irregular  frequencies  are  those  for 
which  there  exist  piecewise  continuous  nontrivial  solutions  of 


r(p)  v(p)  ds^  -  0 


for  p  e  CQ  (47) 


where  we  have  replaced  a(p)u(p)  by  v(p).  Adopting  the  notation  of  Kleinman  and 
9 

Roach  we  Introduce  the  boundary  integral  operators  on  C 


Kv  s»J*  v(q)  ds^ 


for  p  e  CQ  (48) 


and 
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v(q)  &|E*al  da<j 


for  ptfio  (49) 


so  that  Equation  (47)  may  be  written 


(I+K*)  v  -  0 


for  p  e  C 


In  order  to  proceed  we  recall  the  jump  properties  of  single  and  double  layer 
potentials:  if  u  is  continuous  on  Cq  (actually  L^(Cq)  suffices),  then 


i±  J  u(q)  Y(p,q)  dsq  -  +  u(p)  +J  u(q)  3^ 


£kt  3nn 

p^-e  p 


(51a) 


lim±  f  u(q)  ds  -  +  u(p)  +  f  u(q)  ds£ 

tr*-C~  •*  a  **  -  a 


(52a) 


where  T77±  means  that  p  approaches  C  from  D. .  Since  C  may  have  corners.  Equations 

O  i  o 

o 

(51a)  and  (52a)  hold  almost  everywhere  on  Cq;  see  Lemma  4,  Appendix  B).  The  fact 
that  Y  Is  given  by  Equation  (34)  rather  then  merely  being  the  fundamental  solution 
Yq  -  -  1/ (2tt | p— q | ) •  does  not  alter  these  properties,  since  y  -  YQ  is  regular  in  D+, 

D  ,  and  CQ.  That  Cq  is  not  closed  is  of  no  consequence  provided  Cq  is  piecewise  Lia- 
punov  and  satisfies  a  cone  condition.  If  we  further  restrict  Cq  so  that  it  is  made 
up  of  segments  lying  on  Liapunov  surfaces  of  Holder  index  1  (|&q*(p~q)|<A|p-q|^), 
then  Leona  6,  Appendix  B  assures  that  if  one  of  the  derivatives 


b\  •<*>  d,<  or  77-f u<”  *£5  d*< 


exists  for  p  e  CQ,  then  the  other  does  also  and  they  are  equal.  Ve  will  henceforth 
assume  that  this  slightly  stronger  smoothness  requirement  is  satisfied,  which  still 
allows  for  corners  and  edges. 

It  will  be  useful  to  add  the  jump  conditions  for  layers  defined  on  C  .  If  u  is 

1  ^  * 
continuous  on  (again  L  (C  )  would  suffice),  then,  since  q  *  q  (see  Equation  (37)), 

the  singularity  at  q  is  doubled  in  strength  and 


Y(p.q)  ■ - - - +  0(!n|p-q|) 

"Ip-ql 

hence 


lim_ 

p*C 

w 


J  u(q)  Y(p,q>  d»q 

c 

w 


-  2u(p)  +  j  u(q)  d8 

C  q 

w 


(51b) 


1  »<«> 

P^w  C  q 

w 


2u(p)  + 


u(q) 


3Y(p,q) 

3”, 


(52k) 


where  p  approaches  from  D_,  since  we  have  not,  as  yet,  defined  Y(p*q)  outside  the 

slab  -h  <_  y  (y  )  _<  0.  The  normal  on  C  is  taken  to  point  In  the  positive  y 
—  p  q  —  w 

direction. 

Another  way  to  arrive  at  boundary  Integral  equations  for  the  floating  body 
problem  is  to  assume  that  a  solution  is  given  in  the  form  of  a  single  layer  with  un¬ 
known  continuous  density  on  Cq. 


♦  (p)  *• 


u(q)  Y(P»q>  <*»„ 


for  p  e  D,  (53) 


Then,  with  Leona  1  in  Appendix  B  we  aee  that  $  is  continuous  to  the  boundary  (we  may 
define  u  -  0  on  3D+\Cq).  Also,  <j>  satisfies:  Laplace's  equation.  Equation  (lb);  the 
free  surface  condition.  Equation  (lc);  the  boundary  condition  on  Cfi,  Equation  (Id); 
and  the  radiation  condition.  Equation  (If).  Moreover,  with  Lemma  9  of  Appendix  B  we 
can  show  that  <J>  t  H^qc  (D+).  Finally,  to  satisfy  the  boundary  condition  on  Cq  we 
take  the  normal  derivative,  use  the  jump  relation.  Equation  (51a),  and  arrive  at  the 
integral  equation 


u(p)  + 


j> 


(,)  d  . 
P 


for  p  a.e.  on  C  (54) 
o 


or,  with  Equation  (48), 


(I+K)u 


for  p  e  Cq  (55) 


Following  Kleinman  and  Roach,  we  call  those  values  of  k  for  which  (I+K)u  ■  0  has 

nontrivial  continuous  solutions  irregular  or  characteristic  values  of  -R  and  values 

of  k  for  which  (I+K*)u  *  0  has  nontrivial  solutions,  we  call  irregular  or  character- 

istic  values  of  (-K*).  Actually,  Kleinman  and  Roach  (see  also  Angell  and  Kleinman 

2 

for  clarification)  treated  solutions  in  L  (C  )  for  smooth  closed  C  in  which  case  it 

°  2  ° 

is  known  that  solutions  of  (I+K)u  »  0,  which  lie  in  L  (CQ),  also  are  continuous. 

In  the  present  case  we  will  consider  continuous  or  piecewise  continuous  solu¬ 
tions  in  order  that  normal  derivatives  of  double  layer  distributions  satisfy  Lemma  6 
(Appendix  B).  If  K  is  compact,  then  these  irregular  values  are  the  same.  However, 

K  is  not  compact  in  the  space  of  continuous  functions  on  CQ  if  the  surface  is  not 
smooth.  Nevertheless,  some  of  the  features  of  compact  operators  are  retained.  To 
see  this  we  first  define  the  adjoint  floating  body  problem  as  the  following  interior 
problem:  Find  <p  e  C  (D_)  with  interior  boundary  flow  or  in  H  (D_)  such  that 


for  p  e  D  (56) 


3*  vT  ■  n 

37  '  *  0 


for  p  £  Cw  (57) 


♦  -  f 


for  p  e  Cq  (58) 


where  f  is  a  given  function. 

Further,  those  values  of  k  for  which  a  nontrivial  <t>  exists,  such  that  Equations 
(56)- (58)  hold  with  f  *  0,  will  be  called  eigenvalues  of  the  adjoint  floating  body 
problem. 

We  now  prove: 

Theorem  1  -  The  eigenvalues  of  the  adjoint  floating  body  problem  are  real. 

Proof  -  Assume  u  is  an  eigenfunction.  Then  Green's  theorem  implies 


f  (u  ^  - 5  Is) ds  - 0 

r 


_  c 

o|j  w 


(59) 


and,  applying  the  boundary  conditions,  we  get 


I 


(k-k) 


I  u| ^  ds 


'21J 


Im  k 


I  u I ^  ds 


(60) 


hence,  either  u  *  0  or  1m  k  *  0.  However,  if  u  ■  0  on  Cw>  then  since  u  also  vanishes 
on  Cq,  u  is  an  eigenfunction  of  the  interior  Dirichlet  problem  and  it  is  known  that 
there  are  no  nontrivial  solutions  of  this  problem  for  the  Laplacian.  Hence,  Im  k  *  0 
which  proves  the  theorem.  Next  we  prove: 

Theorem  2:  If  k  is  an  eigenvalue  of  the  adjoint  floating  body  problem,  then  k  is  a 
characteristic  value  of  -K. 

Proof:  Assume  k  is  such  an  eigenvalue  so  that  there  exists  a  nontrivial  function  v_ 
defined  in  D_  which  satisfies  Equations  (56),  (57),  and  (58)  with  f  •  0.  With 
Green's  theorem.  Equation  (40),  we  have 


for  p  e  D_ 

(61) 

for  p  e  D 


» 


21 


where  the  boundary  conditions  on  Cq,  as  well  as  the  fact  that  both  y  and  v_  satisfy 
the  free  surface  condition  on  Cv  have  been  used.  Now,  taking  the  normal  derivative 
and  using  the  jump  conditon  of  Equation  (51a)  we  obtain 


2 


9v 

3n 


9v_ 

7T 

q 


ds 

q 


for  p  e  C  a.e.  (62) 
o 


or,  with  Equation  (48), 


3v 

(I+K)  -  0 


for  p  e  C  (63) 
o 


Since  both  v  and  3v  /9n  cannot  both  be  zero  on  C  because  v  is  assumed  a  non- 
—  —  o  — 

trivial  solution  of  the  Laplacian  in  D_,  it  follows  that  k  is  a  characteristic  value 
of  -K.  Furthermore,  we  can  also  show: 

Theorem  3:  If  La  k  >  0  and  k  is  a  characteristic  value  of  -K,  then  k  is  an  eigen¬ 
value  of  the  adjoint  floating  body  problem. 

Proof:  Assume  (I+K)<J>  *  0  has  a  nontrivial  solution  4>.  Then  define 

v+  :■  j  <Kq)  y(p,q)  ds^  for  p  G  D+  (64) 

o 

The  jump  condition.  Equation  (51a) ,  then  shows  that 

3v+ 

m  (I+K)  <J>  “  0  (65) 

and  hence,  the  Uniqueness  Theorem  guarantees  that  v+  -  0  for  p  e  D+. 


Now  let 


Since  the  single  layer  is  continuous  across  Cq  (Lemma  1,  Appendix  B),  It  follows 
that  v_  *  0  on  Cq  and,  because  of  the  properties  of  the  Green's  function,  either  v 
is  an  eigenfunction  of  the  adjoint  floating  body  problem  or  v_  *  0  for  p  e  D  .  The 
latter  possibility  is  ruled  out  since,  if  v_  =  0  then  3v_/3n-  «  0  and,  with  Equations 
(51a)  and  (48),  we  would  have 


3v_ 

■  -4>  +  Rt>  *  0  for  p  e  CQ  (67) 


With  Equation  (65)  this  would  imply  $  ■  0  contrary  to  the  assumption,  thus  proving 
the  theorem.  Theorems  2  and  3  establish  the  equivalence  between  eigenvalues  of  the 
adjoint  floating  body  problem  and  characteristic  values  of  -K.  We  can  also  obtain 
some  relations,  though  not  as  complete,  for  characteristic  values  of  -K*  as  follows. 
Theorem  4:  If  Im  k  >  0  and  k  is  a  characteristic  value  of  -K*,  then  k  is  an  eigen¬ 
value  of  the  adjoint  floating  body  problem,  hence  also  a  characteristic  value  of 
(-K). 

Proof :  Let  ^  be  a  nontrivial  solution  of  (I+K*)  $  *  0  and  p  e  CQ.  Such  a  <f>  exists 
since  k  is  assumed  to  be  a  characteristic  value  of  -K*. 

Define: 

v_  :*  J*  <Kq)  ds^  for  p  e  D 

Co 


Then,  using  the  jump  condition.  Equation  (52a),  we  see  that  v_  ■  (I+K*)  <)>  *  0. 

Hence,  because  of  the  properties  of  y,  either  v_  is  an  eigenfunction  of  the  adjoint 
interior  problem  or  v_  =  0.  The  latter  case  is  ruled  out  by  the  following  argument. 
If  v_  *  0  for  p  e  D_,  then  3v_/3n”  *  0  on  Cq. 

Now  define 


!"  j  4><q)  dsq 


for  p  e  D+ 


and  use  the  continuity  of  the  normal  derivative  of  the  double  layer  (Lemma  6, 
Appendix  B)  to  deduce  that 
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0 


for  peC 
o 


Then  v+  *  0  for  p  e  D+  from  the  Uniqueness  Theorem  and  in  particular,  from  Equation 
(52a), 


(-I+K*)  <|>  -  0 


for  p  e  C 


Since  it  is  also  true  that  (I+K*)  $  ■  0,  it  follows  that  4>  “  0,  contrary  to  assump¬ 
tion,  thus  establishing  that  k  is  an  eigenvalue  of  the  adjoint  floating  body  problem. 
Application  of  Theorem  2  completes  the  proof. 


UNIQUELY  SOLVABLE  INTBGRAL  EQUATIONS 

We  saw  in  the  previous  section  that  the  integral  equation  formulation  of  the 
problem  led  to  equations  which  were  not  uniquely  solvable  at  irregular  frequencies. 
Here  we  present  two  methods  for  modifying  these  equations  to  regain  unique  solvabili¬ 
ty.  One  method  involves  modifying  the  domain  of  the  operator,  whereas  the  second 
involves  modifying  the  operator  on  the  same  domain. 

Method  1;  This  method  involves  the  use  of  Equation  (38)  not  only  on  Cq,  but  in  D_ 

as  well.  In  scattering  theory  this  is  sometimes  called  the  extended  boundary  con- 

21  22  23 

dltion  method  (see,  e.g.,  Kupradze,  Schenck,  and  Waterman  and  has  been  used  in 

slightly  different  form  in  the  present  context  by  Chang  and  Pien. ^  First  we  prove: 

Theorem  5:  If  u  is  continuous  and  bounded  on  C  and 
-  o 


J. 


uM£aSLd9 

dn  q 


for  p  e  D_lJCw  (68) 


then  u 
Proof:  Let 


0  on  C  . 
o 


"  f  u  dsc 

•J  a 


for  p  e  D_yCw 
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Then, 


-  0 


for  p  e  D  C 
r  -u  w 


and,  using  Equation  (52a) 


and 


llm  <p 

p-C~ 

o 


0 


“(p)  +  J  u(q) 
C 

o 


3®,  i 


for  p  a.e.  on  C 
r  o 

(69) 


3$ 

-  ■  0  for  p  a.e.  on  C 

3n”  ° 


Now  define 


♦+ 


ZjlEiSl 

dn 

q 


for  p  e  D+ 


With  Lemma  6,  Appendix  B,  on  the  continuity  of  normal  derivatives  of  the  double 
layer  with  continuous  density,  we  have 


for  p  e  C 

o 

hence,  4>+  is  a  solution  of  the  homogeneous  floating  body  problem  and  by  the  Unique¬ 
ness  Theorem,  $+  ■  0  for  p  €  D+.  Again,  with  Equation  (52a)  we  have 

-u(p)  +  J*  u(q)  48^  *  o  for  p  e  CQ  a.e.  (70) 

C  q 

o 
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Subtracting  Equation  (70)  from  Equation  (69)  establishes  the  theorem.  This  means 
that 


f  ♦(ohjjml  * 

r  9 


J  Y(P.q)  V(q)  d8q 

C 

o 


for  p  e  D  C 
r  -u  w 


(71) 


has,  at  most,  one  solution,  but  the  equation  still  presents  problems  since  the  domain 

and  the  range  do  not  coincide  and  the  range  is  all  of  D  C  .  However,  we  can  im- 

-U  w  * 

prove  things  somewhat  by  applying  Equation  (38)  only  on  the  boundary  of  D_  and 
establish 

Theorem  6:  If  u  is  continuous  and  bounded  on  C 


u(p)  +  f  u(q)  ds  ■  0 

*_  n  ^ 


for  peC'  (72) 


and 


u(q)  ds  »  0 

”9  q 


for  peC'  (73) 


then  u  *  0  for  p  e  C  ,.C  .  We  use  C'  and  C'  to  denote  the  smooth  parts  of  C  and  C  . 

oU  w  o  w  o  w 

Proof:  Assume  u  is  a  solution  of  Equations  (72)  and  (73).  Define 


*_(P) 


u(q)  8T<P»1?  ds 
q 


for  p  e  D_  (74) 


Then,  since  u  is  assumed  continuous,  with  Equations  (52a),  (72),  and  (73)  we  see  that 
<t>_  —  0  for  p  e  Co(JCw.  Hence,  -  0  for  p  e  D_  since  there  are  no  elgenfunctione 
of  the  interior  Dirichlet  problem  for  Laplace's  equation.  Then  the  previous  theorem 


shows  that  u  "  0  on  Cq.  This  theorem  shows  that  even  though  k  may  be  an  Irregular 
frequency,  which  means  that  Equation  (72)  may  have  nontrivial  solutions,  by  adding 
the  additional  requirement  of  Equation  (73),  then  uniqueness  is  assured  at  all  fre¬ 
quencies.  Thus,  the  inhomogeneous  integral  equations  in  the  form 


and 


<J>(p)  +J  $(q)  dsq 

C  q 


V(q)  Y(p,q)  dsq 


for  p  e  C' 
o 


(75) 


dn  q 

q 


v(q)  y (p,q)  dsq 


for  p  e  C'  (76) 


have,  at  most,  one  continuous  bounded  solution,  <P (p)  for  p  e  C  .  He  defer  a  proof 

o 

of  existence  of  a  solution  of  these  equations.  These  equations  still  suffer  from 

the  drawback  that  the  unknown  function  $  is  defined  on  CQ,  but  the  equations  must  be 

satisfied  on  C  .  ,C  .  A  more  usual  equation  of  the  second  kind  with  a  kernel  whose 
oU  w 

domain  and  range  are  the  same  may  be  obtained  through  the  use  of  the  following. 
Theorem  7:  If  7T(p)  is  continuous  and  bounded  on  C^,  and 

2  "u(p)  +  I  TT(q)  dsq  -  0  for  p  e  Cw  (77) 

JC  q 

w 


then  u(p)  ■  0.  Here  fi  points  from  C  in  the  positive  y-direction. 

q  * 

Proof;  Assume  7  satisfies  Equation  (77)  and  let 


$ 


for  p  e 


Then,  because  of  the  properties  of  y»  v<  *««  that  <p  satisfies  Laplace's  equation  in 
the  slab,  -h  <  y  <  0;  the  free  surface  condition.  Equation  (lc),  on  Cf;  the 
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homogeneous  Neumann  condition.  Equation  (Id),  on  CB;  the  radiation  condition.  Equa¬ 
tion  (If);  and,  because  of  Equation  (77)  and  the  Jump  condition,  Equation  (52b), 

$  -  0  for  p  £  Cw>  _ 

Now  apply  the  divergence  theorem  to  ♦  over  D+(^0yD_  using  the  above  re¬ 
lations  obtaining 


(78) 


Now  we  follow  precisely  the  part  of  the  uniqueness  proof  following  Equation  (16), 
since  <j>  defined  above  is  of  the  form  of  Equation  (17),  and  conclude  that 

for  for  P  e 


Moreover,  since  y  satisfies  the  free  surface  condition  on  y  -  0  we  also  have 


*  -  0 


u(q)  Y(p,q)  dsq 


for  p  E  D_  (79) 


or 


Juvq)  Y(P»q>  <*sq  “  0 

C 

w 


for  p  £  D_  (80) 


Now  take  the  normal  derivative  form  D_  using  the  Jump  condition.  Equation  (51b), 
obtaining 


-2  ”u(p)  +  f  uXq)  dsq  -  0  for  p  e  Cw  (81) 

•*C  P 

w 

But  with  the  representation  of  the  Green’s  function.  Equation  (34),  we  see  that, 
for  p  and  q  e  Cw  (=^  yp  -  yq  -  0) 
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3y_  3 

3n~  “  3_ 

q  p 


St  "■'£ 


k  (k^-k^) 

—  v —  (k  R)  sixth  k  h  cosh  k  h 

i  on  n  n 


(82 


n-0  hkn-hk+k 


Thus,  Equation  (81)  can  be  written 


-2  u(p)  +  j  'u(q)  ds  "  0 

•l  a  " 


for  p  e  C'  (83 
w 


which,  with  Equation  (77),  implies  that  TT(p)  -  0  for  p  £  C^,  thus  establishing  the 
theorem.  This  theorem  shows  that  if  Equation  (77)  holds,  then  necessarily 

for  p  e  CQ  (8< 


Now  combine  Equations  (75),  (76),  (83),  and  (84)  and  define 


$(p) 


<j»(p)  for  p  e  Cq 


(85 


uT(p)  for  p  e  Cw 


to  obtain  simply 


*(p)  +  J  B(q)  *(q)  dsq  -  J  V(q)  Y(P,q)  ds 


C  C 
ollw 


q 

for  p  e 


(86 


C'  C' 

oUw 


where  8(q) 


1  for  q  £  C 

o 

1/2  for  qeCy 


2' 


Equation  (86)  Is  of  the  desired  fora  and  if  it  has  a  unique  solution,  then  the  re¬ 
striction  of  *  to  C  is  the  solution  of  Equations  (75)  and  (76).  It  remains  then 
o 

to  verify  that  Equation  (86)  has,  at  most,  one  solution  which  is  the  content  of: 

Theorem  8:  If  $  is  piecewise  continuous  on  C_,,C.  and 
-  o  ou  w 


VP)  +J  e(q)  Vq)  d8q  •  0  for  P  e  CoUCw  (87) 


C  C 
oUw 


then  4>  *  0 

o 


Proof:  Assume  ♦  satisfies  Equation  (87)  and  define 
-  o 


v  :-  f  B(q)  $  (q)  h&lll  ds  for  p  e  D  (88) 

-  J  o  3n  q 


C  C 
oU  w 


Then,  using  the  jump  conditions.  Equations  (52a)  and  (52b)  and  Equation  (87),  we 
find  that  v  *  0  for  p  e  3D_.  Hence,  v_  •  0  for  p  e  D_  since  the  interior  Dirichlet 
problem  has  no  eigenfunctions.  Thus, 


3v_ 

JTm  0 


for  p  e  Cq  (89) 


Now  define 


:-  B(q)  *  (q)  ds  for  p  e  D+  (90) 

J  n  “ 


c  c 

oUw 


Since  v+  is  a  double  layer  with  piecewise  continuous  density,  we  may  invoke 
Appendix  B,  to  see  that 


3v+  3v_ 


for  p  e  CQ  (91) 


■.  •:  v  -i  ;s 
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I 


But  then  v+  is  a  solution  of  the  homogeneous  floating  body  problem  and  the  Uniqueness 
Theorem  then  implies  that 


v+(p)  -  0  for  p  e  D+  (92) 

Taking  limiting  values  as  p  -*■  Cq  we  have,  with  Equation  (52a), 

v_  -  0  for  p  e  D_  =&  4>o  +  j  &(q)  $>o(q)  ds  *  0  (93) 

c  c  ^ 

oU  w  for  p  c  CQ 

ds  -  0  (94) 

q 

for  p  e  Cq 

from  which  we  deduce  that  4*0  for  p  e  C  .  Using  this  result  in  Equation  (87)  we 

o  o 

have 

2  *o(p)  +  *o(q)  ds  -  0  for  p  G  Cw  (95) 

C  « 

w 


and 


0  for  p  e  D .  =o  -♦  + 
r  +  o 


/ 


B(q) 


C  C 
oU  w 


But  then  the  unique  solvability  proven  in  Theorem  7  implies  $Q(p)  "  0  for  p  e  . 
Hence,  $Q(p)  "  0  for  p  e  Cq(J  C^,  thus  establishing  the  theorem. 

To  summarize  the  results  of  this  section:  If  (p)  is  a  solution  of  the  floating 
body  problem  Equations  (la)-(lf),  then  on  Co,  $(p)  -  $(p)  for  p  e  CQ  where  ♦ 
satisfies  the  equation  ~ 

$(p)  +  J  B (q)  Hq)  dsq  -  J  V(q)  y(p,q)  dsq  (96) 

CoUCw  Co  for  p  g  C'  C' 

oUv 

and  this  equation  has,  at  most,  one  solution. 
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An  alternate  form  of  the  equation  may  be  derived  whan  V(p)  is  the  normal  deriva 
tive  of  an  Interior  potential;  i.e.. 


V(p)  -  - 


for  p  e  Cq  (97) 


For  example,  in  heaving  motion  V  «  -n  •  y,  hence  <J>  *  y  +  C.  When  Equation  (97) 
holds,  we  may  use  the  Green's  identities.  Equations  (38)  and  (40),  to  write 


I  J  [<Kq)-t?(q)l  -  7  (p*q>  4-  t<t(q)4(q)]  )  ds 

JC  '  q  q  J 

o 

-  2  ?(p)  -  a(p)  (?(p)-H.(p)]  +  J  J^Y(p,q)  - 


T,ns  ilXEiSl" 


dsq  (98) 


Now  define 


4>t(p)  4»  (p)  +  (p) 


and  use  the  boundary  conditions  on  Cq  and  to  obtain  the  representation 


o(p)  <f>t(p)  +J  ^(q)  dsq  -  2  ?(p)  +  J  Y(p,q)  l 

r  9  n  L 


§r  -  k*(q)l  d8c 

q  J  ^ 


(100) 


The  representation  in  Equation  (100)  is  valid  for  all  points  in  the  slab  -h  <_  y  <_  0. 

The  major  advantage  of  this  form  lies  in  the  case  of  heaving  motion  where  we 
may  cause  the  integral  over  the  waterplane  to  vanish  by  choosing  the  constant  in 
$  appropriately,  namely 


Then,  Equations  (75)  and  (76)  became 


♦  (P)  +J  *t(q>  d«*q  -  2  (y-  i) 


j 

•»_  q 


3l 


ds 


q  k 


for  p  e  C'  (102) 


for  p  e  C£  (103) 


whereas,  in  place  of  Equation  (96)  we  have. 


*(P>  +  f  8 (q)  *(q)  5y|S^-  *B  «  2  (y-  £) 

J  Q 


(104) 


c  ,,c 
oU  w 


for  P  e  C^c; 


where 

$(p)  -  4>t(p)  for  P  e  c0 

Method  2:  A  second  method  of  obtaining  an  integral  equation  with,  at  most,  one  solu¬ 
tion  which  does  not  involve  extending  the  domain  of  the  integral  operator  is  pat- 
terned  after  the  method  used  in  acoustic  scattering  problems  by  Burton  and  Miller, 
Kleinman  and  Roach,9  and  Angell  and  Kleinman.10  Again  we  start  with  the  representa¬ 
tion,  Equation  (38);  evaluate  it  on  Cq;  take  the  normal  derivative  from  D+,  using 
Equations  (51a)  and  (52a);  and  obtain  the  pair  of  equations 

<J>(p)  +  J  <Kq)  d8q  -  J  v<9>  f(p.q)  d8q  for  P  e  co  (105) 

C  q  C 

o  O 

and 

J-  f  <t>(q)  d,  .  f  V(q)  ds  -  V(p)  (106) 

3np  J  3nq  q  J  %  q 

0  °  for  p  E  C' 
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f 


f 


Since  V  is  assumed  to  be  in  L^(Cq)  then  Lemma  5,  Appendix  B,  ensures  that  the  normal 
derivative  of  the  right-hand  side  of  Equation  (105)  exists  in  and  since  3$/3n  is 
also  assumed  to  exist  (3$/3n  is  equal  to  V),  then  two  of  the  three  terms  in  Equation 
(38)  have  normal  derivatives  from  D+,  hence  the  third  term  also  must  have  a  normal 
derivative.  In  addition,  since  <J>  is  assumed  continuous.  Lemma  6,  Appendix  B,  en¬ 
sures  that  the  normal  derivative  of  the  double  layer. 


3 

3n 

P 


I 


<P(q) 


p>q) 

q 


exists  almost  everywhere  on  and  is  the  same  whether  we  approach  Cq  from  D+  or  D  . 
Hence  Equation  (106)  is  obtained  even  though  Cq  may  have  corners.  Denoting  this 
normal  derivative  by  Dn,  i.e., 

V!’5rf  d*q  <107> 

p  c  q 

o 

we  may  rewrite  Equations  (105)  and  (106)  as 

(I+K*)<J>  -  J  V(q)  y(p,q)  dsq  for  p  e  C'  (108) 

C 

o 

D  <j>  -  KV  -  V  for  p  e  C'  (109) 

n  o 

Alternately,  if  V(p)  is  the  normal  derivative  of  an  interior  potential.  Equation 
(97),  then  we  may  differentiate  the  representation  of  Equation  (100)  in  the  normal 
direction,  obtaining 


3 

3n 


|  Vq)  ~T'^2'<L)'  d8q  "  "2  V<P> 


.  f  3Y(p,q> 

J  3nn 


(110) 


and  the  pair  of  boundary  integral  equations  become 


(I+K*)  A  -  2  <}>(p)  + 


J  Y(p.q)  -W(q)J 


|| - k£(q)J  dsq  for  peC'  (111) 


Dn*t  = 


-2  V(p)  + 


I  ^  (I-  -*->) 

«  p  \  q  / 


|| - k$(q))  ds^  for  p  e  C'  (112) 

c  p  ^  q 


w 


As  with  Equation  (100),  we  remark  that  this  form  is  convenient  in  the  case  of 
heaving  motion  with  <J>  =  y  -  1/k  in  which  case  the  waterplane  integrals  vanish.  It 
is  the  pair  of  equations,  either  Equations  (108)  and  (109)  or  Equations  (111)  and 
(112),  which  will  be  shown  to  be  uniquely  solvable.  First  we  observe: 

Theorem  9:  If  k  is  not  an  eigenvalue  of  the  adjoint  floating  body  problem  and  is 
a  piecewise  continuous  solution  (in  the  closure  of  Cq)  of 

( I+K* )<f>  =  J  V (q)  y(p,q)  d8q  for  p  e  C'  (113) 

C 

o 


then 


D  <f>  -  KV  -  V 
n 


for  peC'  (114) 


Proof:  Let  4>  satisfy  Equation  (108).  Then  define 


/_  :=J  4>(q)  dsq  -  J  V (q)  y(p,q) 


dsq  for  p  e  D_  (115) 


and  use  the  jump  condition.  Equation  (52a),  to  see  that  v_  ■  0  on  Cq  and,  in  fact, 
that  either  v  »  0  for  p  e  D_  or  v  is  an  eigenfunction  of  the  adjoint  floating  body 
problem.  But,  since  k  is  not  an  eigenvalue  of  this  problem,  v_  ■  0  for  p  e  D_  and 
hence  has  a  normal  derivative  from  D__.  The  term 

J  V (q)  Y(p,q)  dsq 

C 

o 
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also  has  a  normal  derivative,  Equation  (51a),  hence  the  normal  derivative  of  the 
double  layer  exists  and  the  theorem  follows. 

Thus,  if  k  is  not  an  eigenvalue  of  the  adjoint  floating  body  problem  it  is  suf¬ 
ficient  to  solve  Equation  (108).  A  similar  argument  shows  that  under  the  same 
assumption  on  k,  it  suffices  to  solve  Equation  (111). 

It  remains  to  show  that  regardless  of  the  value  of  k,  the  pair  of  Equations 
(108)  and  (109)  or  Equations  (111)  and  (112)  are  uniquely  solvable.  In  this  section 
we  concern  ourselves  only  with  uniqueness  and  prove 

Theorem  10:  If  Im  k  ^  0,  the  only  piecewise  continuous  solution  of  the  pair 


(I+K*)<J>  -  0 

(116) 

o 

i 

(117) 

is  <J>  =  0. 

Proof :  Assume  4>  satisfies  Equations  (116)  and  (117)  and  construct  the  function 

v+  :»  I  <p( q)  ds  for  p  e  D+  (118) 

C  9  q 

o 

Since  Equation  (117)  is  satisfied,  it  follows  that 

for  p  e  C 
r  o 

Hence,  v+  *  0  for  p  £  D+,  otherwise  it  would  be  a  nontrivial  solution  of  the  homo¬ 
geneous  floating  body  problem.  How  take  the  limit  as  p  approaches  Cq  from  D+  using 
the  jump  condition.  Equation  (52a),  to  obtain 

lim+  v+(p)  -  -4>(p)  +  j  4>(q)  dsq  “  0  (119) 

P^o  JC  q 

o 

This  last  equation  may  be  written  as 


-*■ .  ■' 
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(-I+K*)<j>  -  0 


for  p  e  Cq  (120) 

which,  when  used  with  Equation  (116)  implies  that  $  *  0. 

As  a  consequence  of  this  result,  it  follows  that  Equations  (108)  and  (109)  (or 
Equations  (111)  and  (112))  have,  at  most,  one  solution.  Furthermore,  Theorem  9  shows 
that  when  k  is  not  an  eigenvalue  of  the  adjoint' floating  body  problem  and  <J>  is  a 
solution  of  Equation  (108)  it  is  also  a  solution  of  Equation  (109).  Similarly,  if 
<J>  satisfies  Equation  (111)  it  also  satisfies  Equation  (112). 

The  question  of  existence  has  not  yet  been  settled,  however,  we  conclude  this 
section  with  one  more  uniqueness  result  which  is  advantageous  from  a  numerical  view¬ 
point  in  that  it  involves  a  single  integral  equation  in  contrast  to  a  pair.  This  is 

O 

patterned  after  a  result  for  the  Helmholtz  equation  by  Burton  and  Miller.  The  idea 
is  to  combine  Equations  (108)  and  (109)  by  multiplying  one  of  the  equations  by  a 
suitable  constant  (or  function)  r)  and  adding  it  to  the  other,  thus  obtaining  a 
single  equation 

(I+K*+nDn)4>  -  J  V (q)  y(p,q)  dsq  +  n  (KV-V)  for  p  e  Co  (121) 

C 

o 

Remarkably  this  equation  has,  at  most,  one  solution  If  Im  n  J1  0.  We  state  this  as 
Theorem  11;  If  Im  n  t  0,  Im  k  *  0,  and  iji  is  a  piecewise  continuous  solution  of 

(I+K*+nD  )<f>  =  0  for  p  e  C  (122) 

n  o 


then  <|>  ■  0. 

Proof:  Assume  $  satisfies  Equation  (122)  and  define 


v 


3y(p.q) 

q 


ds 

q 


for  p  e  D_  (123) 


With  the  jump  condition.  Equation  (52a),  and  Equation  (123)  we  have 
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v  -  (14**)$  -  -n  D  $ 

—  tl 


for  p  £  C(  (124) 


and 


8v-  1 

_T  ,  Dn$  -  -  i  (14**)$  for  p  e  Co  (125) 


Since  y  satisfies  the  free  surface  condition  on  C^,  v_  will  also  satisfy  this  con¬ 
dition.  Applying  Green's  theorem  on  D  we  get 


0 


J  (v- 

C  ,C  ' 
oU  w 


3v 

9n 


(126) 


where  Equations  (124)  and  (125)  have  been  used.  The  integral  over  vanishes  since 
Im  k  *  0,  hence  v  and  v  satisfy  the  same  free  surface  condition  on  C^.  Since 
ri  -  ri  ^  0  on  Cfll  Equation  (126)  implies  that 

Dn$  -  0  for  p  e  Cq  (127) 

which,  with  Equation  (122),  also  Implies  that 

(14**)$  -  0  for  p  e  Cq  (128) 


Now  define 


/  «i>  il5£iai  d», 
C  ** 


for  p  e  D+  (129) 


With  Equation  (127)  and  the  continuity  of  the  normal  derivative  of  the  double  layer 


1 


Thus,  v+  ■  0  for  p  E  D+  otherwise  It  would  violate  the  Uniqueness  Theorem.  Taking 
the  limit  as  p  approaches  Co  from  D+,  and  using  Equation  (52a),  we  find 


lim  v  -  (-1+K*)(J>  -  0 

p+C 
r  o 


(131) 


Equations  (128)  and  (131)  then  imply  that  <p  -  0  as  required  and  the  proof  is 
complete. 

In  summary,  if  4>  (p>  is  a  solution  of  the  floating  body  problem.  Equations  (la)- 
(lf)  then,  on  Co,  <J>(p)  satisfies  the  equation 


*(P>  *<■»  +  "  TST  !  »<’>  ■% 

C  q  P  c  q 

o  o 

-  J  V (q)  jjy (p,q)+H  -Y^n,A)-  dsq  -  r)V(q) 


(132) 


and  this  equation  has,  at  most,  one  solution  if  Im  n  £  0- 

Alternately,  if  V(p)  -  -3T/3n  for  p  e  Cq  where  4»  is  a  solution  of  Laplace's 
equation  in  D_,  and  4>t(p)  -  T(p)  +  <J>(p)»  then  4>t  satisfies  the  equation 

♦t<P)  +  J  Vq)  d8q  +  n  3ST  j  Vq)  d8q 

J  o  o  _  q 


-  2  ?(p)  +  J  ^Y(p»q)+n  d8q  "  2qV<P) 


(133) 


for  p  e  C 


and  this  equation  has,  at  most,  one  solution.  In  the  case  of  heaving  motion,  when 
V(p)  -  -ny  *  -8  •  9*  this  becomes 
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<J»t(p)  +  f  *t(q>  d8q  +  *  f  Vq)  d8q  "  2  (y~  £)  +  2n  n 

**  fl  n  *  11  ^ 


y 

(134) 


Having  established  the  uniqueness  results  our  task  is  not  complete  until  we 
show  that  the  integral  Equations  (96)  or  (132)  actually  do  have  a  solution,  however, 
this  existence  proof  is  not  included  here. 
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APPENDIX  A 

A  TRANSCENDENTAL  EQUATION 

We  list  here  some  properties  of  the  roots  of  the  equation  defining  the  wave 
frequencies 


k  sinh  k  h  =  k  cosh  k  h 
n  n  n 


(A.l) 


Property  1:  If  kn  is  a  root,  then  -kn  is  a  root. 

Property  2:  If  la  k  -  0  and  Re  k  ^  0,  then  there  exist  exactly  two  distinct  real 
roots  and  an  infinite  number  of  purely  imaginary  roots.  Thia  was  observed  by  Frit* 
John. 

Property  3:  If  Im  k  ^  0,  then  Re  k^  4  0.  This  may  be  seen  by  assuming  that  the 
statement  is  not  true  i.e..  Re  kn  -  0  or  that  kn  *  ixn,  in  which  case.  Equation  (A.l) 
becomes 


-x  sin  x  h  *  k  cos  x  h 
n  n  n 


(A.  2) 


Since  the  left-hand  side  is  real  and  Im  k  ^  0,  we  find  cos  x  h  *  0  in  which 

n 

case,  from  Equation  (A. 2),  x  sin  x  h  *  0  and  these  two  equations  cannot  be  true 

n  n 

s imultaneously . 

Properties  1  through  3,  considered  together,  establish 
Property  4:  If  k  f  0,  there  exists  at  least  one  root  of  Equation  (A.l)  with  a 
positive  real  part. 

Property  5:  If  k  is  bounded,  then  Re  kR  is  bounded.  This  may  be  seen  as  follows. 

Assume  k  is  a  root  of  Equation  (A.l),  such  that  |k  | >2 1 k| .  Then,  solving  Equation 
k  h  n  “ 

(A.l)  for  e  n  we  see  that 


2k  h  k  +k  -2k  h  k  -k 
n  n  .  n  n 

6  "  k~Tk  and  e  "  k+k 

n  n 


(A.  3) 


from  which  it  follows  that 


41 


(A. 4) 


2h|RekJ 

e 


|kj+|k| 

|kj-|k| 


Since  the  right-hand  side  is  monotonically  decreasing  in  (kj  for 


2h  Rek 


2 1  k  1  k  1 
2 1  k|  —  |  k| 


or 


Recalling  that  this  was  found  under  the  assumption  that  | >  2| 
either  Equation  (A. 6)  holds  or 


IkJ  <  2 jkj 


which  implies 

|Rekn]  <  2|k| 


thus. 


|Re^nl  I 


max 


|2|k|,  ^  *n3 


which  establishes  Property 
Equation  (A. 4), 


Property  6: 


lim  Rek 


n’ 


This  helps  establish 


5.  Not  only  is  Rekn  bounded, 

-  0 


but,  as 


|kni  2|k|  we  have 
(A.  5) 

(A.  6) 

c|  we  see  that 


(A.  7) 

can  be  seen  from 


Property  7:  There  exists  a  root  of  Equation  (A.l)  with  largest  positive  real  part. 
We  designate  the  root  as  kQ. 

This  follows  since  any  bounded  sequence  of  roots  is  either  finite,  hence  has  a 
member  with  largest  real  part,  or  converges  to  a  bounded  root  and  again  there  is  one 
with  largest  real  part.  The  bound  may  be  chosen,  with  Property  6,  so  that  the  roots 
with  magnitude  larger  than  the  bound  have  real  parts  smaller  than  Rek  .  Finally  we 
add 

Property  8:  If  Rekn  >  0  and  la  k  >  0,  then  Im  kQ  •>  0.  This  can  be  seen  by  letting 
kn  -  xn  +  iyn  with  xr  >  0  and  rewriting  Equation  (A.l)  as 


k  -  k 


slnh  k  h 
_ n 

n  cosh  k  h 
n 


x  sinh  2x„  h-y^  sin  2y  h+i(y  sinh  2x  h+x  sin  2y  h) 
n _ n  •'n _ ^n  Jn _ n  n _ -^n 

2 | cosh  (xn+iyn)hj? 


(A.  8) 


But,  Im  k  >  0  implies  that  y  sinh  2x  h  +  x  sin  2y  h  >  0.  Recall  that  x  >  0 
—  n  n  n  n  —  n 

and  assume  that  yn  <  0.  Dividing  by  2h  xn  yn  yields 


sinh  2x  h  sin  2y  h 

“  +  -  —  -  <  0 


2x  h 
n 


2ynh  - 


hence. 


sinh  2x  h  sin  2y  h 

n  <  _  —  in  < 

2x  h  -  2y  h  - 

n  'n 

which  is  impossible  with  xn  4  0.  Hence,  yn  >  0  establishing  the  desired  result. 


sin  2yn  h 


2y  h 
'n 


<  1 
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APPENDIX  B 


SOME  PROPERTIES  OF  SINGLE  AND  DOUBLE  LAYER  POTENTIALS  ON 
NONSMOOTH  BOUNDARIES 


n 

If  C  is  a  simply  connected  piecewise  Liapunov  surface,  that  is,  C  -  U  C 

i«l  iU 

n  ••  13 

U  3  C. ,  where  each  C.  lies  on  a  Liapunov  surface  o.  (e.g.,  Gunter  ),  then  many  of 

i*l 

the  familiar  properties  of  single  and  double  layer  potentials  on  smooth  boundaries 
are  still  valid  as  shown  by  Wendland . ^  We  list  some  of  these  useful  properties  here 
with  some  indication  of  how  they  are  obtained.  First  we  define  the  symbols: 


V 


V 


V 


c 


c 


y(q)  Y0(p,q)  dsq 


3YQ(p.q) 

y(q)  __5_ - d 

q 

3YQCp.q) 

v(q)  -  ds 

p 


for  p  e  IR3  (B.l) 


for  p  e  F?  (B.2) 


(B.  3) 

for  p  e  for  some  i 


where 


Y„(P,q)  5-  ~ 


2tt  t  p-q  f 


(B.4) 


We  understand  to  denote  an  open  set  on  C  so  that  corner  and  edge  points  are  not 
Included  in  C^,  but  are  confined  to  SC^.  Note  that 

D  v  *  K*y  for  p  e  C 


where 


9Y  <P.q) 

v(q>  — —  dsq 


for  p  c  C  (B.S) 
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While  y  is  defined  only  on  C,  we  extend  by  0  onto  a^,  that  is,  define 

y1(p)  :■  M(p)  for  p  E  0^3^ 

(B.6) 

0  for  p  e  ot\(CiJiCi) 

Also,  if  is  the  Bolder  index  of  the  surface  a^,  let  a  -  min  {ot^}  so  that  a  may 
serve  as  a  common  index  for  all  a^.  1 

Lemma  1:  v  e  l*(C)^Sn  t  C°,a  (|^)=>  S  y  e  LP(C)  for  p  >  1. 

°  10  ° 

Proof:  Angell  and  Kleinman  show  that  the  result  holds  if  C  is  Liapunov  with  index 
a.  But, 


V 


z\ 

i-1  C, 


y(q)  Y0(p,q)  <uq 


si 

i-1  o. 


(q)  Y0(p»q)  dsq 


(B.  7) 


i-1 


where 


Soi  yi 


^(q)  Y0(p.q)  dsq 


(B.8) 


But 


y  e  L° (C)  ^yjCL*  (ot) 

hence 

Sol  «1  E  c°'“  (K3’ 


for  each  i  as  does  a  finite  linear  combination,  thus  establishing  the  result 


Lemma  2:  y  e  L*  (C)  Kqu  e  LP(C)  for  p  ^  1, 

Proof:  As  before,  we  rewrite  the  operator  as  a  sum  of  operators 


^(p.q) 


v-i’J  pi<<>)  — % 

1-1  c±  q 


(B.9) 


Now,  for  any  p  e  R“ 


3Y_(p.q) 

<q>  d\ 


,  fin* (p-q) I 

1 2w  l!pH  »  |  -~9 — rj— d® 

1  L  (C Id-o  3 


L  (C) 


(B.10) 


-c  i  p-q| 

ci 


Rewrite  the  integral  on  the  right  as 


/ 


|n  •(p-q)l 

5 —  da 

3  q 


I  p-q  I 


■[ 


cinVp) 


IV  <p-q)  I 

— 9 - =—  ds  + 

I p-q 1 3  q 


I 


cA(cinVp» 


|nQ*(p-q)| 

- 3  d8q 


l  p-q  l 


(B.ll) 


where  (p)  is  a  ball  of  radius  6  with  the  center  at  p,  and  we  choose  <S  to  be  less 
than  half  the  Liapunov  radius  of  the  Liapunov  surface  of  which  is  a  part. 
Since 


l p-q I  >  $ 


(B.12) 

for  q  e  Ci\JCin  Bfi(p)) 


\<cmVp» 


lfiq*(p-q)l 


%<J 

ciVcinV*» 


ds  H(C  ) 

|p-q|2l~ 


(B.'3) 


where  M(C1)  is  the  surface  area  of  C^.  To  evaluate  the  first  Integral  on  the  right 
of  Equation  (B.ll)  we  consider  two  cases.  If  p  e  C^c  o^,  there  exist  constants  A 


and  a  such  that 


|fiq*(p-q>l  1  A  |p-ql 


(B.14) 


hence 


f  |n  *(p-q)| 

— 9 - =—  ds 

J  ,  v  I p-q I  3  q 

CinVP>  1Pql 


r  ds 

<  A  — a~ 

-  Jcln»5W  lp-’> 


5  „2ir 


<u/  <pj 


dpJ  d<|>  <  A'6a 

00  P  (B.15) 


for  some  A'  where  we  used  the  standard  method  of  projection  onto  the  plane  tangent 
to  at  p.  Here  A'  is  independent  of  p  and  6.  If  p  ^  then,  denoting  by  p*  the 


normal  projection  of  p  onto  o^,  we  have 


I"  *(p-q>! 


|n  • (p-p+p-q) ] 


cmVp> 


— -  ds 

|3  q 


CinVp) 


-= -  ds 

1 3  q 


cin®6(p)  lp~qi  C.  _  B 


I (p-q) 1 


~d8q 


m  5 


, (B.16) 


Note  that  |p-‘p|  <  6  even  if  p"  i  Ci  and  that  B^(p)  c  B2g(p).  Using  this  fact  and 
Equation  (B.14)  which  holds  with  'p  replacing  p. 


f  |n  *(p-q)j 

__a — __  d8 

L'*<» q 


Cin  26tp; 


f  Ip-fll1^  ds  (B.17) 

W>'-'  ’ 


Now  we  use  the  inequality  (e.g.,  Gunter  ) 


p-q|  (Ip-pI  +|p-q|2)1/ 


(B.18) 


to  show  that 


■  :4  Z*-4?*  * 


Jzm 


Employing  Equations  (B. 18)  and  (B.19)  in  Equation  (B.17),  we  obtain 


1 


|8  ‘(p^-q) 

— * -  ds 


clnB6<>>) 


lp-q 


A  Ud 

3  q 


1  2 1 P-?l  J 


ds 


^nB2S(P>(lp-?|2+l?-q|2)3/2 


ds 


+  84  J  _  Tr-fcr 

I™' 


25  2tt 


<  2  j  p 


-? i  f  dp  r 

J  Jr 


26  2ir 


d<j> 


0  ■'o 


-.2^3/2  +  8Af  “4 

x'fi  Jn  M 


(|p-'p|  +p) 


'0  v0 


<  - 


<IP-?IV>1/2 


25 


+  ~  (26 )“ 


(B.20 


p”0 


from  which  we  infer  that  there  are  constant  A^  and  independent  of  p  and  6  such 
that 


1 

cinV'>) 


1p„- (p-q)  1 


ds  <  A.  +  B. 
q  -  i  i 


a 

O 


(B.21 


Although  Equation  (B.21)  was  derived  assuming  that  p  $  C^,  we  see,  with  Equation 
(B. 15),  that  Equation  (B.21)  remains  valid  even  if  p  e  C^.  Employing  the  estimates 
of  Equations  (B.13)  and  (B.21)  in  Equation  (B.10),  we  have 


Setting 


1_  max 
"  2tt  i 


M(Ci>  R  .a 
~2 ~  +Ai+Bi  6 


(B.23) 


we  have  finally,  from  Equation  (B.9),  that 


i  *  i 

K u 


n 

i«i 


J 


9y  (p»<i) 

-  d8q 


<  nM  y 


L  (C) 


(B.24) 


oo  ft  9t  n 

Thus,  £  L  (C)  =>  KqM  is  bounded,  hence,  KqM  c  1/ (C)  for  p  j>  1,  which  concludes  the 
proof  of  Lemma  2.  By  Fubini's  Theorem  we  then  obtain 
Lemma  3:*  y  e  LP(C)=>  KQy  e  L^CC)  for  p  ^  1. 

Next  we  have 

Lemma  4:  If  y  e  L^(C),  then 


lim  -  S  u  =  +  y(p)  +  K  u 
?  on  o  o 

p-C+  p 


for  p  a.e.  on  C 
(B.25) 


and 


lim_  DQy  -  +  y(p)  +  KQy 

P-C+ 


for  p  a.e.  on  C 
(B.26) 


*We  remark  that  K  y  need  not  be  bounded  (e.g. ,  Fichera  and  Snelder-Ludovlci, 


24 


25 


26, 


Leis,  and  Craggs,  Mangier  and  Zamir  ). 


Here  the  +  denotes  limits  from  either  side  of  C.  We  choose  +  to  denote  the  portion 

o  3 

of  *  into  which  the  normal  from  C  points,  call  it  D+;  then  D_  -  JR  \p+  C 
Proof:  Since  the  results  are  to  hold  almost  everywhere,  corner  end  edge  points  of 
C  may  be  excluded  and  if  p  is  a  smooth  point  of  C^,  then 


lim  -r —  S  y  *  li® 

-  3n  o^  •  __ 

p-*C+  P  p+C  j-1  P  C 


M(q)  Y(p,q) 


j 


n  '  f 

•Zh 


3Y0(P»q)  a  f 

(,)  -fj -  d.  +  lta  ^  I  p(q)  Y„(p.q>  d.q 

1-1  P  P*C+  'C, 


**  ^  r  3y  (P*q)  3  f 

Z  — d8, +  lta,  ta-  J  viM  r<,(p'q)  d‘< 

J  p  p-C+  P  o. 


j-1  C 


j 


(B.27) 


where  I'  means  that  the  i  term  is  omitted. 


Since  is  Liapunov,  the  validity  of  the  jump  relations  (Gunter  )  shows  that 

"  r  3y  (p»q)  _  f  ,  a  ,  ,  , 

L^V-2,  J  p<">  -fe - dVi(p)  +J  «4h)  55-  % 

^  P  i  n  P  fT 


lim 

p^C 


j-1  c. 


f  8Yn(P»q)  _  f  N  3  ,  .  , 

-  y\  v<q)  — Is — d8  +y(p)+J  M(q)  airYotp^)  d*q 


j-l  C 


J 


3Y_(p»q) 

y(q)  “1^ -  d*q  +  W(P) 


(B.28) 


A  similar  decomposition  establishes  the  Jump  condition  for  the  double  layer. 
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*  *.« 


As  a  consequence  of  Lemmas  2,  3,  and  4,  we  have 
Lemma  5:  If  y  e  L*(C),  then 

llm-  1HT  Soy  e  l1(C)  (B.29) 

p-C+  P 

i 

i 

I  and 

ilm_  DQy  e  LP(C)  for  p  >  1  (B.30) 

p-*C+ 

This  follows  immediately  since  y  £  L* (C)  implies  that  Lemma  4  holds  and  also 

Lemmas  2  and  3,  hence  y  and  KQy  are  in  LP(C)  and  Kqu  e  Ll(C).  Also  of  interest  is 

the  normal  derivative  of  the  double  layer  for  which  we  may  state 

igggg-ii  If  y  £  C0^)  for  every  i,  a±  is  Liapunov  with  index  1,  (which  means 

|n  •(p-q)|<A|p-q|  ),  and  if  one  of  the  derivatives  3D  y/3n+  or  3D  y/3n~  exists  at  a 
“  o  o 

point  p  £  C,  then  the  other  exists  and  they  are  equal. 

£ro.Pl.:.  He  obviously  exclude  corner  and  edge  points  since  there  is  no  normal  at  such 
points.  If  p  is  an  interior  point  of  C  ,  then  the  proof  in  the  smooth  case  (Gunter 

13  1 

p  297  et  eq.  )  may  be  repeated  without  change  since  the  crucial  step  Involves  treat¬ 
ment  of  a  small  patch  around  p  which,  in  this  case  as  in  the  smooth  case,  lies  on  an 
appropriately  smooth  Liapunov  surface.  The  same  argument  shows  that  the  result  holds 
for  piecewise  continuous  densities  even  if  the  discontinuities  occur  at  smooth 
points,  provided  all  points  of  discontinuity  are  excluded. 

Finally,  we  list  some  properties  of  single  and  double  layers  considered  as 
operators  mapping  functions  defined  on  C,  not  to  functions  on  C,  but  to  functions 
defined  on  the  interior  of  C  which  we  denote  as  D_.  From  Leans  1  we  immediately  have 
Lemma  7:  y  e  L* (C)  *  Soy  £  LP(D_)  for  p  >  1. 

Also,  since  D  y  £  C  (D  )  and  lim  Due  L°°(C) 

p-*C 

(Lemmas  2  and  4),  for  y  e  L* (C),  we  have 
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Lena  8:  y  e  L  (C)=*  D  u  e  LP(D  )  for  o  >  1. 

— — — —  o  —  — 

More  difficult  to  establish  are  properties  of  derivatives  of  potentials.  However, 
we  do  have 

oo  0 

Lena  9:  y  e  L  (C)  =»  7  SQy  e  L  (D_) 

Proof :  For  p  e  D_,  we  have 

l7Vli4rj  l“(Dl  I"  (».3i> 

C  |p-ql  L  (C)  J  |p-q| 


.therefore. 


J  hVi2dt  i-i-riiPii2.  j 

i  P  (2it)2  L  (C)  J  p-q  2  J 


p-q|  £  lP' 


i— ^ iiuii2.  fff 

(2ir)Z  L  (C>  J 


(2n)2  iTinJ  J  J  I - TTT - ^d8od»o  <B*32> 

(2ir)  L  (C)-',  JQ  ^  |  p-q | 2 1  p-q^ | 2  •»  «»i 


.2* _ <  f 

1 2 j _ _  ,2  ~  J 


D_  iP-'ll  l^ll  Bd(0)  lp|  fp^^ll' 


(B.33) 


where  d  is  twice  the  diameter  of  D_.  Explicitly 


.d  .it  .  2tt 


j  .  U  *  II  4» 

f  - r2 - 2  if  drf  d6f  *”7 - SfLi - 

JD  |p-q|2|p-q1|2  rZ-f|q-qi!  -2p*  (q^q) 


f  I  q-qi  *  r27r  ff  »  l2 

<  J  dr  d0  dd>  sin  e  ^  ^  Pn(ooa  y)  ♦ 

0  J0  0  Ln-0  'q“qll 


<».»> 

(coot.) 


\ 


where  we  have  employed  the  generating  function  for  Lengendre  polynomials  and 


p*<q,-q) 

cos  y  .  - ± -  (B.  35) 

Ip^qj-q)! 


Using  the  orthgonality  of  Legendre  functions  leads  to 


<  8tt  1  m  12tt 

”  n«o  (2n+l)2  )q~qxl 

Substituting  this  estimate  in  Equation  (B.32)  we  obtain 

i 


(B.  36) 
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f  Hvl2*pifllull’  f  I  ^ 


«n  c  i«i 


'i|2.,  f 

L  (C)J 


ifllullV  J  SG1  dsq 


CB.  37) 


But  Lemma  1  ensures  us  that  the  single  layer  with  density  1  is  bounded 
Integral  over  C,  hence 


as  Is  its 


I  |7v|2*pi!iw^(c)  "vnil(c) 


(B. 38) 


which  completes  the  proof  of  the  Leona. 

It  should  be  noted  that  no  claim  is  made  about  the  sharpness  of  the  results  of 
this  Appendix.  Indeed  more  complete  and  precise  estimates  of  the  mapping  properties 
of  single  and  double  layers  may  become  available  (see,  for  example,  Fabes,  Jodelt, 
and  Riviere^). 
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